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Abstract 

Two-dimensional hydromagnetic flow of an incompressible Jeffrey nanofluid over an exponentially stretching surface is 
examined in the present article. Heat and mass transfer analysis is performed in the presence of thermal radiation, viscous 
dissipation, and Brownian motion and thermophoresis effects. Mathematical modelling of considered flow problem is 
developed under boundary layer and Rosseland's approximations. The governing nonlinear partial differential equations are 
converted into ordinary differential equations via transformations. Solution expressions of velocity, temperature and 
concentration are presented in the series forms. Impacts of physical parameters on the dimensionless temperature and 
concentration are shown and discussed. Skin-friction coefficients are analyzed numerically. A comparison in a limiting sense 
is provided to validate the present series solutions. 
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Introduction 

Boundary layer flow over a stretched surface is a subject of 
abundant studies now a days because of its existence in various 
engineering and industrial processes like cooling of metallic sheets 
in a cooling bath, annealing and thinning of copper wires, 
aerodynamic extrusion of plastic and rubber sheets, drawing of 
plastic films and sheets, glass fiber and paper production etc. It is 
worthmentioning to point out here that the stretching velocity is 
not linear necessarily in all the cases. The stretching velocity may 
be nonlinear or exponential. For example in annealing and 
thinning of copper wires, the desired quality product depends on 
the continuous stretching of surface with exponential dependence 
velocity distribution [16]. The magnetohydrodynamic flow is 
quite interesting in engineering and industrial processes. El 
Koumy et al. [7] investigated peristaltic flow of Maxwell fluid 
through a porous medium in the presence of Hall effects. Hall 
currents and heat transfer analysis in peristaltic flow is performed 
by Abo-Eldahab et al. [8] . Effects of magnetic field and porous 
space in peristaltic flow of Maxwell fluid are examined by 
Mekheimer et al. [9]. Shehzad et al. [10] discussed the Joule 
heating and thermophoresis effects in MHD flow of Jeffrey fluid 
induced by a stretching sheet. An induced magnetic field and slip 
effects in peristaltic flow through a porous medium are described 
by Mekheimer et al. [1 1]. Recently, Hayat et al. [12] performed a 



study to analyze the effects of Hall current and Ohmic heating in 
peristaltic flow of non-Newtonian fluid. 

A working fluid is involved in many engineering and industrial 
applications that are in flowing or stagnant state. This working 
fluid is used to transfer energy/heat from one position to the other. 
The adequate heat transfer performance has been a major issue 
for a long period. This issue can be resolved by using a new 
working fluid that has better thermal performance than the 
ordinary base liquids. Recently nanofluid is the best candidate to 
take place of working fluid. Nanofluid is a fluid in which the 
nanoparticles are submerged. The size of these nanoparticles is 1— 
1 00 nm. The thermal conductivity of the nanofluids is higher than 
that of base fluids. Further, the novel properties of Brownian 
motion and thermophoresis of such fluids make them potentially 
useful. Nanoparticles are used to enhance the thermal character- 
istics of ordinary base fluids such as water, ethylene glycol or oil 
[13]. In addition the magneto-nanofluid is a unique material that 
has both liquid and magnetic properties. Such nanofluid has 
superficial role in construction of loud speakers, blood analysis and 
cancer therapy. Buongiorno [14] provided a mathematical model 
of nanofluid which has the characteristics of thermophoresis and 
Brownian motion. Later on, Makinde and Aziz [15] investigated 
the boundary layer flow of viscous nanofluid with convective 
thermal boundary condition. Here, the flow is induced due to 
linear stretching of surface. Closed form solutions of MHD 
nanofluid flow with heat and mass transfer analysis in the presence 
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of slip condition were developed by Turkyilmazoglu [16]. Ibrahim 
and Makinde [17] analyzed the effects of thermal and concentra- 
tion stratifications in boundary layer nanofluid flow by a vertical 
surface. Second law of thermodynamics in MHD steady flow of 
nanofluid over a rotating disk was discussed by Rashidi et al. [18]. 
Moradi et al. [19] presented the series solutions of Jeffrey Hamel 
flow by considering the different types of nanoparticles. Unsteady 
natural convection flow of nanofluid with heat and mass transfer 
over a vertical plate was examined by Turkyilmazoglu and Pop 
[20] . Slip and Joule heating effects in MHD peristaltic flow of 
nanofluid under thermal diffusion and diffusion thermo effects was 
studied by Hayat et al. [21]. Sheikholeslami et al. [22] carried out 
an analysis to discuss the effects of an applied magnetic field in 
rotating flow of nanofluid with heat transfer. 

The aim of present article is to study the flow analysis of Jeffrey 
fluid [23-28] in the presence of thermophoresis, Brownian motion, 
thermal radiation and viscous dissipation effects. The flow caused 
is by an exponentially stretching sheet. Jeffrey fluid has ability to 
exhibit the properties of ratio of stress relaxation to retardation 
and retardation. Mathematical formulation is performed under 
boundary layer and Rosseland's assumptions. Homotopy analysis 
method (HAM) [29—35] is utilized for solution expressions of 
dimensionless velocity, temperature and concentration. Temper- 
ature and concentration fields are shown and discussed for the 
different values of arising parameters. Skin-friction coefficients are 
computed and analyzed. Comparison of local Nusselt number in a 
limiting sense is tabulated and analyzed. 

Mathematical Model 

We consider the two-dimensional hydromagnetic flow of Jeffrey 
nanofluid over an exponentially stretching sheet. Heat and mass 
transfer effects are taken into account. An applied magnetic field of 
strength Bq is encountered normal to the flow direction. The 
magnetic Reynolds number is chosen small and the Joule heating 
effects are absent. Further the induced magnetic field is smaller in 
comparison to the applied magnetic field and is negligible. In 
addition, viscous dissipation effects are also taken into account. 
The two-dimensional boundary layer equations of an incompress- 
ible Jeffrey nanofluid with heat and mass transfer are given below: 



du dv q 

dx dy 



(i) 




tl f , kg, % 

Figure 1. h — curves for functions f(rj) , 9(rj) and (/)(rj) at 21th 
order of approximations when X \ = 0.3, ft = 0.2, M= 0.5, Pr 
0.7. Le= 1.0, Nt= 0.2= Nb, Ec= 0.3= Rd, A= 0.1and£= 0.2. 
doi:1 0.1 371/journal.pone.01 0371 9.g001 
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The boundary conditions for the considered flow analysis are 



u= U w (x)= £/ 0 exp(y), v = 0, T=T w = T ao + To 
(Bx 



exp 



Ax\ 



C=C w = C oo + C 0 exp( — J at j = 0, 



(5) 



Table 1. Convergence of homotopy solution for different order of approximations when k \ = 0.3, 
Pr = 0.7= Le, N t = N b = 0.2, Ec= 0.3= Rd, A= 0.1, B= 0.2 and h f = - 0.5= h e = . 



P = 0.2, M= 0.5, 



Order of 



approximation 


- /" ( 0) 


- 6 ' ( 0) 


- r ( 0) 


1 


1 .2900 


0.5316 


0.5572 


5 


1.4310 


0.2972 


0.3806 


11 


1.4350 


0.2433 


0.3146 


20 


1.4350 


0.2210 


0.2692 


30 


1.4350 


0.2124 


0.2435 


45 


1.4350 


0.2076 


0.2225 


50 


1.4350 


0.2076 


0.2193 


60 


1.4350 


0.2076 


0.2193 
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Figure 2. Variation in dimensionless temperature (?(//) vs r\ for 
different values of h when p = 0.2, M= 0.5, Pr = 0.7, 
Le= 1.0, Nt= 0.2= Nb, Ec = 0.3= Rd, A= 0.1 and B = 0.2. 
doi:1 0.1 371 /journal, pone.01 0371 9.g002 



w->0, T^ T f , C^-Cqo, when j>->oo 



(6) 



where u and v are the velocity components in the x — and 
y — directions, v the kinematic viscosity, X\ the ratio of 
relaxation to retardation times, A2 the relaxation time, p / the 
density of fluid, o the Steffan-Boltzman constant, a the thermal 



diffusivity, T 



( PC) p 



the ratio of nanoparticle heat capacity and 



( pc) f 

the base fluid heat capacity, Db the Brownian diffusion coefficient, 
Dj the thermophore tic diffusion coefficient, q r the radiative heat 
flux, Too and Coo are the ambient fluid temperature and 
concentration far away from the sheet and A, B, Tq, Co are the 
constants. 

By employing Rosseland's approximation, Eq. (3) has the form 
dT dT 



dx dy 

16a*r^ d 2 T 



du 



dCdT 



dy dy 
VI2 du d 



D^fdT 

Too \dy 
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du du 



The dimensionless variables are defined as 



u=U 0 exp (y) f'(rj), v = - y ^ exp (f (rj) + rjf'(rj)), 

^ = ^ exp S)' r=r - +r ° exp (f) m - 



(8) 




Figure 4. Variation in dimensionless temperature 0(ij) vs r\ for 
different values of Pr when Ai=0.3, P = 0.2, M= 0.5, 
Le = 1.0, Nt= 0.2= M>, £c = 0.3= Rd, A= 0.1 and B= 0.2. 
doi:1 0.1 371/journal.pone.01 0371 9.g004 

The equations of linear momentum, energy, concentration and 
their corresponding boundary conditions in dimensionless form 
can be written as 



/'" + (1 + h)ff" ~ 2(1 + h)f' 2 + P( If" 2 ~ \ff" +ff") 
-M(l+2i)/' = 0, 



(9) 



(1 + AO Pr (f0' - 4f0) + (1 + AO Pr M>0'f 



- (1 + AO Pr M0' 2 + Pr £ c f/" 2 + J/"(3/7" -ff") ) = 0, (10) 



(/)" + Pr W - Bf'$) + (TV,/ TV* ) 0" = 0, (11) 
f = 0,f = 1,0' = !,(/)' = 1 atf/ = 0, 

/'-►0, 0->O,^-»O as 77^00, (12) 



where 



/ 



is the Deborah number, 



(tBqU w . . v 

M = is the magnetic parameter, Pr = — is the Prandtl 



p f c 



number, Ec = 



CpT 0 \ Uq 



4 - A 



is the Eckert number, 





Figure 3. Variation in dimensionless temperature 0(i]) vs y\ for 
different values of M when X x = 0.3, p = 0.2, Pr = 0.7, 
Le= 1.0, Nt= 0.2= Nb, Ec= 0.3= Rd, A= 0.1 and B= 0.2. 
doi:1 0.1 371 /journal. pone.01 0371 9.g003 



Figure 5. Variation in dimensionless temperature 0(i]) vs // for 
different values of Nt when X\ = 03, p = 0.2, M= 0.5, Pr = 0.7, 
Le= 1.0, Nb= 0.2, Ec= 0.3= iW, A= 0.1 and 5= 0.2. 
doi:1 0.1 371/journal. pone.01 0371 9.g005 
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Figure 6. Variation in dimensionless temperature 0(ij) vs y\ for 
different values of Nb when X x = 0.3, p = 0.2, M= 0.5, 
Pr = 0.7, Le= 1.0, Nt= 0.2, £c = 0.3= Rd, A = 0.1 and B = 0.2. 
doi:1 0.1 371 /journal.pone.01 0371 9.g006 



( pc) p D B C 0 exp 



Le = is the Lewis number, Nb 



Bx 
~2f 



( pc) f v 

the Brownian motion parameter, 

Ax 



( pc) D T T 0 Qxp 



Nt-- 



ter, 



21 



( pc) f vTn 



is the thermophoresis parame- 



The skin friction coefficient, the local Nusselt and Sherwood 
numbers are defined below. 



c f 



;,Nu x 



xq w 



,Sh x -- 



xq„ 



PfUlix)' x k(T w -T m y x D B {C W -C X ) 



(13) 



where t w is the shear stress along the stretching surface, is the 
surface heat flux and q m is the surface mass flux. The local skin- 
friction coefficient, local Nusselt and Sherwood numbers in 
dimensionless forms can be written as 



Vz^*C fx = j—r- (f"(0) + pf"(0)\ Nu x /Re l J 2 - 
1 -rM 



(14) 



^/R e y 2 =-^'(0), 

where Re x = U w (x)x/v is the local Reynolds number. 



Figure 8. Variation in dimensionless temperature 9(rf) vs // for 
different values of Rd when X x = 0.3, ft = 0.2, M= 0.5, 
Pr = 0.7, Le= 1.0, Nt= 0.2= Nb, Ec = 0.3, A= 0.1 and 5= 0.2. 
doi:1 0.1 371/journal.pone.01 0371 9.g008 



Homotopy Analysis Solutions 

To proceed the homotopic solutions, the initial guesses and 
auxiliary linear operators are chosen as follows: 

/oO/)=l- exp (-*/), 0 o O/) = Qxp(-f]),(j) 0 (rj)= exp (-*/), (15) 



Uf) =f" -f, UP) = 0" ~ 0, m = r~l (16) 

The above initial guesses and auxiliary linear operators satisfies 
the below mentioned properties 

L(f)(C l + C 2 e 11 + C 3 e- 11 ) = 0, 

L(6)(C 4 e r, + C 5 e- r, ) = 0, (17) 

L(^)(C 6 ^ + C 7 e-*) = 0, 

where C/( i= 1 — 7) are the arbitrary constants. 
The zeroth order problems can be written as 

(1 -</)£(/) [fa; q)-Mn)} =qhfNf\f(m ?)], (18) 
(\-q)L(8)[8( m q)-e 0 (r,)] = ^N 9 [fa; 9 ), fa,?), fa, <?)],(19) 

{\-q)m[Kn\ q)-Mn)] =qme\f(m q),e~(n,q),Kn,q)\,W) 



Figure 7. Variation in dimensionless temperature 0(i]) vs y\ for 
different values of Ec when X x = 0.3, ft = 0.2, M = 0.5, 
Pr = 0.7, Le= 1.0, Nt= 0.2= Nb, Rd= 0.3, A= 0.1 and B= 0.2. 
doi:1 0.1 371 /journal.pone.01 0371 9.g007 



Figure 9. Variation in dimensionless concentration (/>(//) vs r\ 

for different values of X \ when £ = 0.2, M= 0.5, Pr = 0.7, 
Le= 1.0, i\ft = 0.2= Nb, Rd= 0.3= ^= 0.1 and B= 0.2. 
doi:1 0.1 371/journal.pone.01 0371 9.g009 



PLOS ONE | www.plosone.org 



4 



August 2014 | Volume 9 | Issue 8 | e103719 



Radiative Hydromagnetic Flow of Jeffrey Nanofluid 




Figure 10. Variation in dimensionless concentration <;>()/) vs r\ 

for different values of Pr when X\ = 0.3, /? = 0.2, M= 0.5, 
Le = 1.0, Nt= 0.2= #6, JW = 0.3= Ec, A= 0.1 and 5 = 0.2. 
doi:1 0.1 371 /journal.pone.01 0371 9.g01 0 



f(O;q) = Oj(O;q)=\,0(O,q) 
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PrEc/3 (df(r,,q)d 2 f(r,,q) ~ .8 2 f(r,,q) 



drj drj 



(23) 
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Figure 11. Variation in dimensionless concentration (/>(//) vs // 
for different values of Le when Ai = 0.3, £ = 0.2, M = 0.5, 
Pr = 0.7, j\ft = 0.2= M>, Rd= 0.3= ^= O.lo and B= 0.2. 
doi:1 0.1 371/journal.pone.01 0371 9.g01 1 



d 2 (/)(rj, q) 



J>vLeb(n, q ) ^ 



B- 



<K*i, q) 



(24) 



(N t /N b ) 



d 2 0(rj,q) 
drj 2 



where fy, he and h§ are the non-zero auxiliary parameters, 
q e [0, 1] is an embedding parameter and Ny, N# and are 
the nonlinear operators. Putting q= 0 and q= 1 one has 

Am o)=Mr,),d(mO)=eo(ri),kmO)=Mfi) and/fa i)=/fa), 

(25) 

%,1) = %),^,1) = ^). 

When we increase the values of q from 0 to 1 then f(rj, q) , 
0 (rj , q) and ^ ( rj , vary from / 0 (i/), 0 O O7), ^ 0 (ff) to 
f(rj), 6( rj ) and ^ ( 7/ ) . By adopting Taylor series expansion, 
we have [28-30]: 



f(V,q)=Mv)+^fm(ri)<r, 

m = 1 

oo 
m = 1 

CO 

(j>(mg)=<l>o(^+ ^<L(W, 



(26) 

(27) 

(28) 




Figure 12. Variation in dimensionless concentration vs ,7 
for different values of Nt when Ai = 0.3, ft = 0.2, M= 0.5, 
Pr = 0.7,Le= 1.0, Nb= 0.2, iW= 0.3= A= 0.1and£ = 0.2. 
doi:1 0.1 371/journal.pone.01 0371 9.g01 2 



PLOS ONE | www.plosone.org 



5 



August 2014 | Volume 9 | Issue 8 | e103719 



Radiative Hydromagnetic Flow of Jeffrey Nanofluid 




Figure 13. Variation in dimensionless concentration </>(//) vs // 
for different values of Nt when X x = 0.3, p = 0.2, M = 0.5, 
Pr = 0.7, Le= 1.0, Nt= 0.2, Rd = 0.3= Ec, A = 0.1 and B= 0.2. 
doi:1 0.1 371 /journal.pone.01 0371 9.g01 3 




Figure 15. Variation in Nusselt and Sherwood numbers for 
different values of Nb vs Nt when X x = 0.3, p = 0.2, M= 0.5, 
Pr = 0.7, Le= 1.0, Nt= 0.2, Rd = 0.3= Ec, A= 0.1 and 5= 0.2. 
doi:l 0.1 371/journal.pone.01 0371 9.g01 5 



fmin)-- 



1 d m f( m q) 



ml drf 



,0 m (rj)- 



1 d m 6( m q) 



q = 0 



<i>m(n)= 



ml df] m 
1 d m $( m q) 



ml drf 



q = 0 



q = 0 



(29) 



The convergence of above series highly depends upon the 
suitable values of h f, h q and h $ . Considering that h f, h q 
and h $ are selected properly such that (26)-(28) converge at 
q= 1 and then we have 



m = l 



6(rj) = e 0 (ri)+J29 m (ri), 



The general solutions can be written as 



fmin) =PJn) + C\ + C 2 e« + C 3 e~ 



e m (r j ) = e t m (r,)+C 4 e" + C 5 e-', 



(30) 



(31) 



(32) 



(33) 
(34) 



tminWmW+Gn* + 0,6-", (35) 
where f^(rj), 0* m (rj) and ^ m {fj) are the special solutions. 




Figure 14. Variation in Nusselt and Sherwood numbers for 
different values of Nb vs Nt when X x = 0.3, p = 0.2, M= 0.5, 
Pr = 0.7, Le= 1.0, Rd= 0.3= Ec, A = 0.1 and B= 0.2. 
doi:1 0.1 371/journal.pone.01 0371 9.g014 



Convergence Analysis and Discussion 

The auxiliary parameters h /, ho and h $ are encountered 
when homotopy analysis method has been utilized to compute the 
series solutions. These parameters have essential importance for 
adjusting and controlling the convergence of series solutions. The 
appropriate values of these parameters are required for the 
convergent solutions. To obtain the proper values of these 
auxiliary parameters, we drawn the h — curves at 2\ th -order of 
HAM approximations. These h — curves are presented 
in Fig. 1. From this Fig. we examined that the suitable values of 
hf, h e and h& are -0.65< h f <-0A0, -0.70 < h 0 <-030, 
— 65 < hfo < — 0.30. The series converges in the whole region 
of rj when hj = — 0.5= ho = h^ (see Table 1). 

Figs. 2-8 are drawn to examine the variations in dimensionless 
temperature profile 6 ( rj ) for different values of ratio of 
relaxation to retardation times 1\, magnetic parameter M, 
Prandtl number Pr , thermophoresis parameter Nt, Brownian 
motion parameter Nb, Eckert number Ec and radiation 
parameter Rd. Fig. 2 shows that an increase in ratio of relaxation 
to retardation times gives rise to the temperature and thermal 
boundary layer thickness. Minimum temperature and thinnest 
thermal boundary layer thickness is noticed when X \ = 0. An 
increase in X \ implies to an increase in relaxation time and 
decrease in retardation time. This change in relaxation and 
retardation times elucidates the higher temperature and thicker 
thermal boundary layer thickness. Change in temperature profile 
for different values of magnetic parameter is noticed in Fig. 3. 
Here we examined that both temperature and thermal boundary 
layer thickness are enhanced for larger magnetic parameter. It is 
also seen that for M= 0, hydrodynamic flow situation is 
recovered. Magnetic parameter involves the Lorentz force. Higher 
magnetic parameter implies to stronger Lorentz force and lower 
magnetic parameter has weaker Lorentz force. Here the stronger 
Lorentz force leads to an increase in the temperature and thermal 
boundary layer thickness. Fig. 4 depicts the impact of Prandtl 
number on the temperature field. It is observed that lower 
temperature and thinner thermal boundary layer thickness appear 
corresponding to the increasing values of Prandtl number. Prandtl 
number is the ratio of viscous to thermal diffusivities. Larger 
Prandtl number fluids have higher viscous diffusivity and weaker 
thermal diffusivity. Such change in viscous and thermal diffusiv- 
ities leads to a decrease in thermal boundary layer thickness. 
Figs. 5 and 6 illustrate that both temperature and thermal 
boundary layer thickness are enhanced when the larger values of 
thermophoresis and Brownian motion parameters are taken into 
account. From Fig. 7 we noticed that temperature is an increasing 
function of Eckert number. Further Ec= 0 corresponds to the 
analysis when viscous dissipation effects are absent. Effects of 
radiative parameter on the temperature field are analyzed in 
Fig. 8. We observed that an increase in radiative parameter leads 



PLOS ONE | www.plosone.org 



6 



August 2014 | Volume 9 | Issue 8 | e103719 



Table 2. Numerical values of skin friction coefficient 
stretching. 



Radiative Hydromagnetic Flow of Jeffrey Nanofluid 
'2ReCf x for different values of X \ , M and /? in case of exponential and linear 





M 


o 

p 


- y/lReCfy, 


-V2Re Cf x 








Exponential 
stretching case 


Linear 

stretching case 


0.0 


0.5 


0.2 


1.63612 


1.34164 


0.4 






1.38278 


1.13389 


0.7 






1.25485 


1.02899 


1.0 






1.15689 


0.94868 


0.3 


0.0 


0.2 


1.25661 


0.96077 




0.7 




1.50007 


1.25269 




1.0 




1.59258 


1.35873 




1.4 




1.70796 


1.48842 


0.3 


0.5 


0.0 


1.28616 


1.07417 






0.3 


1.50353 


1.22474 






0.5 


1.63159 


1.31559 






0.8 


1.80611 


1.44115 



doi:1 0.1 371/journal.pone.01 0371 9.t002 



to an enhancement in the temperature field. In fact an increase in 
radiative parameter provides more heat to fluid that corresponds 
to higher temperature. 

Change in concentration distribution function (j)(rj) for various 
values of ratio of relaxation to retardation times X\, Prandtl 
number Pr , Lewis number Le, thermophoresis parameter Nt 
and Brownian motion parameter Nb is examined through the 
Figs. 9-13. Fig. 9 clearly indicates that concentration and its 
related boundary layer thickness are increasing functions of ratio 
of relaxation to retardation times. In addition a comparison of the 
Figs. 2 and 9 shows that the impacts of ratio of relaxation to 
retardation times on the temperature and concentration are quite 
reverse. Fig. 10 depicts that concentration boundary layer 
thickness become thinner for larger Prandtl number. Form 
Fig. 1 1 we noticed that an increase in Lewis number shows lower 



concentration field. Lewis number is inversely proportional to the 
Brownian diffusion coefficient. This Brownian diffusion coefficient 
becomes smaller corresponding to the larger values of Lewis 
number. This smaller Brownian diffusion coefficient shows a 
reduction in the concentration and associated boundary layer 
thickness. Fig. 12 presents that the larger values of thermophoresis 
parameter give rise to the concentration and its related boundary 
layer thickness. It is also seen that concentration is at the peak 
when Nt= 1.5 and rj = 2.0. Fig. 13 is drawn to see the variations 
in concentration profile when Nb= 0.1, 0.5, 0.8, 1.2 and 1.5. 
We have seen that the concentration and its related boundary 
layer thickness is decreased by increasing the Brownian motion 
parameter. It is also observed that the concentration profile 
decreases rapidly when Nt= 0.1 but after Nt= 0.5 this reduction 
is very small (see Fig. 13). Figs. 14 and 15 are drawn to examine 



Table 3. Comparison values of — 9 ' ( 0) with Bidin and Nazar [3] for different values of Ec, Pr and Rd when = 0, X\ = 0, 

M= 0, Nt= Nb= 0, A= 1.0 and h e = - 0.5.. 



Rd 


Bidin and Nazar [3] 




Present HAM solutions 






Ec=0.0 






Ec=0.0 








Pr= 7.0 


Pr=2.0 


Pr=3.0 


Pr= 7.0 


Pr=2.0 


Pr=3.0 


0.0 


0.955 


1.471 


1.869 


0.9548 


1.4715 


1.8693 


0.5 


0.677 


1.074 


1.381 


0.6775 


1.0734 


1.3807 


1.0 


0.532 


0.863 


1.121 


0.5337 


0.8627 


1.1213 


Rd 


Ec = 0.2 






Ec=0.2 








0.862 


1.306 


1.688 


0.8624 


1.3055 


1.6384 




0.618 


0.965 


1.229 


0.6183 


0.9653 


1.2286 




0.488 


0.782 


1.007 


0.4889 


0.7818 


1.0067 


Rd 


Ec = 0.9 






Ec=0.9 








0.539 


0.725 


0.830 


0.5387 


0.7248 


0.8302 




0.410 


0.587 


0.696 


0.41 1 1 


0.5870 


0.6963 




0.334 


0.498 


0.606 


0.3355 


0.4984 


0.6055 
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the variations in Nusselt and Sherwood numbers for different 
values Nb vs Nt for both linear and exponential stretching cases. 
Here we noticed that heat and mass transfer rates at wall are 
higher for linear stretching case in comparison to exponential case. 

Table 1 is computed to examine the values of — /"( 0), — 6'(0) 
and -f(0) when ^=0.3, jff = 0.2, M= 0.5, Pr = 0.1 = Le, 
Nt = 0.2= Nb, Ec = 0.3= Rd, A = 0.1, B = 0.2 and 
hf= O.5 = #0=/ty. From this Table, we analyzed that the values 
of —f"(0) converge from 1 1 th order of HAM approximations 
whereas the values of — 6'(0) and — (t>'(0) converge from 45th and 
50th order of HAM deformations, respectively. Numerical values 
of skin-friction coefficient y/2ReC/ x for different values of X\ , M 
and /? of both exponential stretching and linear stretching cases 
are presented in Table 2. Here we examined that the values of 
skin-friction coefficient are smaller for larger X\ while these values 
are increased when larger values of M and fi are encountered. 
Further, it is noticed that the values of skin-friction coefficient are 
larger for exponential stretching case in comparison to the linear 
stretching case. Table 3 provides a comparison study with the 
existing solutions for different values of Rd, Ec and Pr when 
A = B= 1.0, A\ = ft = M= 0. From this Table one can see 
that our present results have an excellent agreement with the 
results of Bidin and Nazar [3] . 

Conclusions 

Radiative hydromagnetic flow of Jeffrey nanofluid over an 
exponentially stretching sheet is studied. Heat and mass transfer 
phenomena are discussed in the presence of viscous dissipation. 
The main observations of present investigation are written below: 
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